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ABSTRACT: We discuss the dependence of running couphngs on the choice of 
regularization method in a general softly-broken = 1 supersymmetric theory. 
Regularization by dimensional reduction respects supersymmetry, but standard di- 
mensional regularization does not. We find expressions for the differences between 
running couplings in the modified minimal subtraction schemes of these two regu- 
larization methods, to one loop order. We also find the two-loop renormalization 
group equations for gaugino masses in both schemes, and discuss the application 
of these results to the Minimal Supersymmetric Standard Model. 



1. Introduction 



Low-energy N = 1 supersymmetry [1] provides an elegant solution to the naturalness 
problem associated with the origin of the electroweak scale. If supersymmetry proves to 
be correct, it is quite possible that the first superpartner to be discovered will be the 
gluino, because of its high production rate at hadron colliders. Since the gluino mass is a 
physical parameter of the supersymmetry-breaking sector, its determination will provide an 
important clue as to the rest of the sparticle spectrum. Measurements of the masses of the 
other superpartner s will then yield tests of the various specific extensions of the Minimal 
Supersymmetric Standard Model (MSSM) which have been proposed over the years. For 
example, in models with soft supersymmetry-breaking terms which are "universal" in the 
sense usually associated with spontaneously broken supergravity [2], the squark masses 
cannot be less than about 0.8 of the gluino mass. The details of the supersymmetric 
spectrum have been studied under a variety of assumptions and constraints on the other 
unknown parameters of the model in [3,4,5] among many others. Further, there are sum 
rules [6] which relate, for example, the masses of the neutralinos and charginos to the 
mass of the gluino, without involving any other unknown input parameters. It is therefore 
useful to work out the predictions for sparticle masses in terms of the running parameters 
of supersymmetric models, including radiative corrections. 

There is a problem of principle, however, in discussing radiative corrections in su- 
persymmetric models. The most popular regularization scheme for discussing radiative 
corrections within the Standard Model, dimensional regularization [7] (DREG), violates 
supersymmetry explicitly because it introduces a mismatch between the numbers of gauge 
boson and gaugino degrees of freedom. The modified scheme known as dimensional reduc- 
tion [8] (DRED) does not violate supersymmetry, and thus maintains the supersymmetric 
Ward identities. In DREG, supersymmetry is violated in the finite parts of one-loop 
graphs, and in the divergent parts of two-loop graphs. This means that the /^-functions 
for a supersymmetric model will be different for the two schemes starting at the two-loop 
level. Also, the running couplings computed in DREG with modified minimal subtraction 
[9] (MS) will differ from those computed in DRED with modified minimal subtraction 
(DR) by virtue of finite one-loop effects. 

In this paper, we give the difference between running couplings in the MS and DR 
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schemes, to lowest non-trivial order, for a general N = 1 supersymmetric model with soft 
breaking. As an application, we then consider the two-loop /^-function for a gaugino mass 
in both schemes. The MS version is obtained by specializing the results of [10,11] for 
a general renormalizable theory to the case of a general N = 1 supersymmetric model. 



The two-loop /3-function in the DR scheme is then derived by simply translating all MS 



couplings into their DR counterparts. We do not enter here into the question of whether a 



fully consistent DR calculation at the two-loop level actually exists; although some results 
have been obtained for supersymmetric theories [12], the use of dimensional reduction in 
non-supersymmetric theories is problematic. However, if such a consistent scheme exists, 
it must reproduce the results obtained here to two-loop order. 

We consider a general N = 1 supersymmetric Yang- Mills model. The chiral superfields 
contain a complex scalar 0j and a two-component fermion which transform as a 
(possibly reducible) representation R of the gauge group G. The superpotential is 

w = (1.1) 

with Y'''^^ = (Y^j]^)*. In addition, the Lagrangian contains soft supersymmetry-breaking 
terms of the form 

-CsB = lh'^''<f>^<|>J<f>k + ^MW + h.c. (1.2) 
where M is the mass of the gaugino A. Strictly speaking, the most general renormalizable 
softly-broken N = 1 supersymmetric model also contains couplings in the superpotential 
with dimensions of (mass) and (mass)'^, and soft supersymmetry-breaking scalar couplings 
with dimensions of (mass)^ and (mass)"^. However, such terms are not relevant to the 
present discussion. 

For simplicity we first give our results for the special case of a simple [or U{1)] gauge 
group. We then explain the modifications required if the gauge group is a direct prod- 
uct in Section 4, where we will also discuss the MSSM. We let = (t^)J denote the 
representation matrices for the gauge group G. Then 

(t^t^)! ^ CiR)S{ 

Tr^(t^t^) = S{R)6^^ 

define the quadratic Casimir invariant C{R) and the Dynkin index S{R) for a represen- 
tation R. For the adjoint representation [of dimension denoted by d{G)], C{G)S^^ = 
jACD jBCD -^{^Y^ f^^^ the structure constants of the group. 

3 



2. Dependence of Running Couplings on the Regularization Method 



The strategy for determining the relationship between couplings in the MS and DR 
schemes is to relate each running parameter to a physical quantity which cannot depend 
on the choice of scheme. 

First, consider the running gaugino mass parameter. By standard techniques, one 
can compute the physical pole mass of the gaugino to one-loop order as a function of the 



running mass parameter in each scheme. In the MS scheme, one finds 



while in the DR scheme, the result is 



In both cases 



1+ 



MS $ 



5C{G) + SC{G)ln^ + ^ 

DR $ 



(2.1) 



(2.2) 



= 25(i?$) J dxxln (^[xm?^ + (1 - x)m^ - x{l - x)M^^^^]/ /j?^ 

and the sum is over all chiral supermultiplets $ = (0, iJj) which couple to the gaugino. 
Comparing (2.1) and (2.2), we obtain (to one-loop order) 

^2 



1 + 



167r2 



C(G) 



(2.3) 



Similarly, we can compute the relationship between the gauge couplings in the two 
schemes. Here we need to make an important distinction between the gauge coupling g 
which appears in the interactions of the vector gauge bosons, and the coupling g which 
occurs in the Yukawa interaction 

jC^ = V2g4>\t^)l{^jX^) + ii. c. 

of the gaugino and the chiral superpartners. Gauge invariance guarantees that g is the 
same everywhere it appears, but only supersymmetry guarantees that g = g. Hence we 
have ^ = ^ in DRED, but we expect ^ g^ in DREG, since the radiative corrections in 
DREG violate supersymmetry. By computing the relevant one-loop graphs in the effective 
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action, and demanding that the physical scattering amphtudes computed in each scheme 
are the same, we find 

r ^2 1 

(2.4) 



9ms ~ 9dr 
and 



Sm=<lm[l + ^IC{G)-C{r)]\ . (2.5) 

Thus there is a different couphng g-^ of the gaugino to each distinct irreducible represen- 
tation r of the chiral supermultiplets. Note that we consistently neglect the distinction 
between g-^^, g-^, and g-^ in the one- loop correction parts of (2.3), (2.4), and (2.5), and in 
the two-loop correction parts of formulas below, since all of the incarnations of the gauge 
coupling are the same to zeroth order. The relations (2.3) and (2.4) have been discussed by 
Yamada [13] (see also [14,15]) for a theory with only gauge vector supermultiplets. Here 
we note that these formulas are not modified by the presence of chiral superfields, which 
is not surprising in view of the fact that DRED and DREG only differ for graphs in which 
there is at least one internal gauge boson line which does not terminate (at either end) in 
a scalar-scalar-gauge boson vertex. 

As a non-trivial consistency check, we consider the two- loop running of both g-^ and 
^MS' by specializing the results of [10,11] for a general renormalizable theory to the case of 
N = 1 supersymmetry. The gauge coupling satisfies 

^|l^-=T^^"+(T6^^" ^'-'^ 

M^) =g^ [SiR) - 3C(G)] 

6(2) =g^ ^-Q[C{G)\^ + 2C{G)S{R) +4S{R)C{R)'j - g^Y'^''Y,jkC{k)/d{G) 

(with all MS couplings on the RHS) as is already known [16]. Here S{R) is the Dynkin 
index summed over all chiral multiplets and S{R)C{R) is the sum of the Dynkin indices 
weighted by the quadratic Casimir invariant. Note that the gauge coupling g^ also satisfies 
eq. (2.6) with gj^ replaced everywhere by gjy^ according to (2.4); this is the well-known 
result that the /^-function for the gauge coupling is scheme-independent through two loops. 

On the other hand, for the coupling of a gaugino to a chiral supermultiplet in the 
irreducible representation r, we find from [11] 



U^) |-10[C(G')]2 + 2C{G)S{R) + hS{R)C{R) - S{R)C{r) - C{G)C{r) + 3[C(r)]2| 



Here g=^S{R) means a sum over all chiral multiplcts of the Dynkin index times the ap- 
propriate g^. It is easy to check that (2.7) is indeed consistent with (2.4), (2.5), and 
(2.6). 

We now consider the Yukawa coupling y^-^^ between a scalar 0j and two chiral fermions 
il^j,'^]^. By again demanding that physical scattering amplitudes computed from the one- 
loop effective action in each scheme are the same, we find 

..^lk f, , g'^ 



MS 



Note that the Y^^'' are totally symmetric in the supersymmetry-respecting DR scheme 
(because of the way they appear in the superpotential) , but not in the MS scheme due to 
the radiative corrections. In A'' = 1 supersymmetry, the scalar quartic interaction 

>Cquartic = -\\/^<P*'<P*^k<Pl (2.9) 
depends entirely on the Yukawa and gauge couplings. In a supersymmetry-respecting 



scheme like DR, this dependence is given simply by 

+ <,2(tA.^*;.+t*Vl) . (2.10) 



\ kl_v. vklm of.\k.\l . \k . 



However, in MS this relation is modified by radiative corrections. Indeed, we find 

again by comparing physical scattering amplitudes derived from the one-loop effective ac- 



tions in the two schemes. This completes the "dictionary" for translating between MS and 



DR couplings at one loop order. A gauge-invariant, supersymmetric fermion mass com- 
ing from a quadratic term in the superpotential will also differ between the two schemes, 
but that difference is trivially derived from (2.8) by taking the scalar with index i to be 
a dummy field with C(r^) = and C{rj) = C{r]^). None of the other couplings (writ- 
ten in component language rather than superfield language) of a general softly-broken 
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N = 1 supersymmetric model differ between MS and DR to this order, because DREG 
and DRED only differ for graphs containing a gauge boson line which does not end in a 
scalar-scalar-gauge boson vertex. 



3. Two-Loop Running of Gaugino Masses 



In this section we will consider the two-loop /^-functions for the gaugino mass computed 
in both DREG and DRED. The results of [11] provide the 2-loop MS /3-function for a 
scalar-fermion-fermion coupling in a general renormalizable theory. From this, it is trivial 
to extract the corresponding result for a fermion mass term, by treating the external scalar 
in eqs. (3.3) and (3.4) of [11] as a dummy field (with no gauge or other Yukawa interactions). 
After further specializing to the case of a gaugino mass in N = 1 supersymmetry, and using 
(2.4) and (2.5) to eliminate g^^ in favor of g^, we obtain: 



The simplest way to find the corresponding equation in DRED is to plug eqs. (2.3) 
and (2.4) into (3.1). This gives 



The method used here avoids potential ambiguities and conceptual problems [12] associated 
with doing two- loop calculations directly in DRED. Wc will return to these issues elsewhere. 
Eqs. (3.1) and (3.2) generalize the expressions given by Yamada[13], who pointed out that 
the one-loop relation = {2M/g)Pg fails to hold at the two- loop level. 




(3.1) 




(3.2) 



4. Direct Product Gauge Groups and the MSSM 



As promised, we now point out the modifications which must be made to the preceding 
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formulas if the gauge group is a product of simple [or U{1)] subgroups Ga- The formulas 
involving radiative corrections for a gauge coupling ga [(2.4)- (2.7)] or gaugino mass Ma 
[(2.1)-(2.3) and (3.1)-(3.2)] for each subgroup G a may be obtained by the following set of 
rules [10]. First, each term which does not involve the quadratic Casimir invariant of a non- 
adjoint representation is diagonal in subgroups; i.e., obtained by g ga, S{R) — > Sa{R), 
and C{G) C{Ga)- For the other terms, one has g'^C{r) ^hd'b^hi'^) (2-5) and 
(2.8); C{k)/d{G) Ca{k)/d{Ga) in (2.6), (2.7), (3.1), and (3.2); and 

g^S{R)G{R) ^ glSa{R) J^dbCbiR) 

b 

in eq. (2.6). [Consistency then determines the replacement rules for the C(r)-dependent 
terms in (2.7).] The correction term in (2.11) becomes a double sum ^h^c9b9c • ■ ■ ^^^^ 
subgroups. Finally, in eqs. (3.1) and (3.2) we need the replacement 

16g^C{R)SiR)M ^ 8glSa{R)J29bCbiR)iMa + M^) . 

b 

In the MSSM, the gauge group is SU{3)c x SU{2)l x U{1)y, with chiral superfields 
Q and L for the S't/(2)i-doublet quarks and leptons, and w, d, e for the 5" t/ (2) i- singlet 
quarks and leptons, and two Higgs doublet chiral superfields Hu and H^. The relevant 
part of the superpotential is 

W = HuQYuU + HdQY^d + H^LYee (4.1) 

where Yu, Y^, Ye are 3x3 Yukawa matrices. The soft supersymmetry- breaking Lagrangian 
includes trilinear scalar couplings 

-C = HuQKu + H^Qh^d + H^LheC + h.c. (4.2) 

where hy j p are again 3x3 matrices in family space, and we use the same symbol for 



scalar components as for the chiral supcrfield. Then in the DR scheme the 2-loop running 
of the three gauge couplings are given by 



^ ^ — 9a 75(1) , 9a 



6=1 x=u,d,e 



(4.3) 



d//^" 167r2 (167r2)2 

Here S^'^ = (33/5, 1, -3) for U{1)y (in a GUT normalization), SU{2)l, and SU{3)c re- 
spectively; Ba^ = (0,4,-18), and 

/ 199/25 27/5 88/5 \ . /26/5 14/5 18/5' 

B'-^I = 9/5 21 24 and C"' = 6 6 2 

V 11/5 9 32 / V 4 4 



The renormalization group equation for the three gaugino mass parameters in DR can then 
be written easily in terms of the same coefficients: 



2^a 



6=1 



x=u,d,e 



(4.4) 



In the MS scheme, the two-loop renormalization group equations for the gauge cou- 



plings are exactly the same as in DR (provided, of course, that one consistently uses 
rather than g^). The MS two-loop renormalization group equations for the gaugino 
mass parameters in the MSSM are given by an equation exactly like (4.4), with the only 
difference that B^a^ must be replaced in (4.4) by B^a^^^ = (0, 16/3, -24). 



5. Conclusion 

In this paper, we have presented analytical expressions for the differences between 
couplings in the MS and DR schemes, to one loop order. The MS scheme, while convenient 
for analyzing the predictions of the standard model, is slightly inconvenient for the MSSM 
because it violates supersymmetry, in the sense that the tree-level supersymmetric relations 
between coupling constants are modified by finite radiative corrections even at one loop. 
We were able to exploit this fact to derive the two-loop /3-function for a gaugino mass in 
DR from knowledge of the corresponding formula in the MS scheme, which in turn followed 
from the results for a general renormalizable theory (using MS) given in [10,11]. The same 
method can in fact be used to derive the renormalization group equations for all of the 
other soft supersymmetry-breaking parameters of the MSSM. We will report on this in a 
future publication. 

Let us close by remarking on the significance of our results for predictions in supergravity- 
inspired models. These are typically obtained by choosing a set of input parameters, in- 
cluding a common gaugino mass, at some very high scale and running the parameters 
down to the electroweak scale. The boundary conditions at the very high scale should pre- 
sumably be applied in a supersymmetry-respecting scheme like DR, and threshold effects 
are simpler [14] in DR. However, it may be more convenient at low energies to work in 
MS. In the MSSM, the two-loop contributions can lower the predicted value of the gluino 
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mass by several percent compared to the one-loop predictions, while having a much smaller 
effect on the other gaugino masses. (In extensions of the MSSM with non-minimal particle 
content, these corrections are potentially much larger.) Another effect which is at least 
as important is the proper treatment of the gluino pole mass, as given by eqs. (2.1) and 
(2.2). One might choose to evaluate the running mass at a scale equal to the running mass 
at that scale (i.e., solve the equation M-^{p) = n or Mjy^(n) = fi, depending on which 
scheme one is working in); however, this will always underestimate the true pole mass of 
the gluino, and often by more than 5 per cent. To roughly estimate the size of the effect, 
suppose the squarks are taken to be degenerate with the gluino and much heavier than all 
quarks; then from (2.2) one finds MpQ^g fa M^{M^)[1 + 90:3/477]. Such considerations 
may have an effect on future efforts to understand the parameters of the MSSM. 
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